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Claims 

A Method for sampling a state space by iteratively generating states xu and their 
weighting factors pi* , wherein the index / is the iteration parameter and the index t 
distinguishes different states xv generated by an iteration /, the method is consisting of 
a first step for selecting an initial sampling distribution function p(x), 
a fifth step for performing an analysis 
and an iteration procedure including 

a second step for generating Nj states xj.t by a numerical sampling algorithm and 
a fourth step for testing at least one criterion to decide whether to continue the iteration 
procedure or to stop the iteration procedure and to go to a fifth step in order to perform the 
analysis using the simulated data, 

characterized in that the iteration procedure further includes 

a third step determining weighting factors pu for states xu generated so far by using 

sampling distribution functions a (*) determined so far and 

a fitting step for determining a sampling distribution function Pj{x) for the next iteration by 
fitting pj(x) to puOfaj) for states xu generated so far, wherein C(xh) is a function, 
respectively a property, of the states xtj . 

Method as claimed in claim 1, wherein the sampling distribution function pj(x) of at least 
one iteration is fitted such that it maximises an objective function preferably defined as a 
function of local comparisons between the sampling distribution function and the product 
puO(xt.t) . 

Method as claimed in claim 1 or 2, wherein the sampling distribution function pj(x) of at 
least one iteration is fitted such that the sampling distribution function pj{x) is large for at 
least one state with a large product puOfas) , and tends to be small for states with a 
small product puCfaj). 

Method as claimed in one of claims 1 to 3, wherein the sampling distribution function 
pj(x) of at least one iteration is a function with at least one constraint, preferably having 
an extreme value in the region of at least one selected state xj^mstr and having vanishing 
values at a distance from the at least one selected state xj^tr , wherein the constraint is 
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preferably a harmonic constraint with at least one constant kcomtr defining the strength of 
the constraint at the Selected state Xj t constr . 

4 

Method as claimed in claim 4, wherein the sampling distribution function pj(x) of at least 
one iteration is the distribution function of the system with constraints or a multicanonical 
distribution function with constraints. 

Method as claimed in claim 1 or 2, wherein the numerical sampling algorithm of at least 
one iteration generates correlated states xu , wherein the sampling distribution function 
pj(x) preferably has a maximum that biases sampling into a sub-region of the state 
space, wherein the sampling distribution function pj{x) is preferably fitted such that the 
maximum is in a region where the product p,\tO{xij) has a maximum, and wherein the 
sampling preferably starts from a state close to the maximum of the sampling distribution 
function pAx). 

Method as claimed in one of claims 4 to 6, wherein the fitting of pj(x) is done by selecting 
states jc fV for which the product p,jO(xij) has extreme values and by using the selected 
states xu to define the region in which pj{x) has extreme values. 

Method as claimed in one of daims 1 to 6, wherein parameters of the sampling distribution 
function pj(x) of at least one iteration are determined by a linear least square fit of the 
logarithm of the un-normalized sampling distribution function pj(x) to the logarithm of the 
product pi,tO(xij) 

Method as claimed in one of claims 1 to 8, wherein the normalization constant of the 
sampling distribution function pj(x) of at least one iteration is estimated from the sampled 
states x fJ and their weighting factors pu. 

Method as claimed in one of claims 1 to 9, wherein at least three iterations are done. 

Method as claimed in one of claims 1 to 10, wherein the function 0(x) is a function of a 
set of at least two functions ®={Oi(x\<Mx\ . ^^(x)} , preferably including at least one 
of the following functions 
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at least one property for which at least one estimate is derived in the analysis of the fifth 
step, , :j 

at least one function thaHs large fbr states that must be sampled to ensure transitions 
between important regions, 

the inverse of the probability distribution function of at least one property of the system, 
and 

the inverse of the probability distribution of the negative logarithm of the distribution 
function of the system . 

A computer software product for sampling a state space by iteratively generating states 
jc,/ and their weighting factors /?/,# , wherein the index / is the iteration parameter and the 
index t distinguishes different states xu generated by an iteration /, said product is 
characterized by a computer-reiadable medium in which program instructions are stored, 
which instructions, when read by a computer, enable the computer to 
select an initial sampling distribution function 
to execute a fifth step for performing an analysis 
and to execute an iteration procedure including 

a second step for generating Nj states xj,t by a numerical sampling algorithm and 
a fourth step for testing at least one criterion to decide whether to continue the iteration 
procedure or to stop the iteration procedure and to go to a fifth step in order to perform the 
analysis using the simulated data, 
characterized in that the iteration further includes 

a third step determining weighting factors for states xu generated so far by using 
sampling distribution functions p { (x) determined so far and 

a fitting step for determining a sampling distribution function pj(x) for the next iteration by 
fitting pj(x) to pmOfaJ) for states xu generated so far, wherein 0(xi.t) is a function, 
respectively a property, of the states xu . 



